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INTRODUCTION 

Detailed design issues associated with y^e^dSension^l, turbulent, vaporizing, 

operation require CFD utility Tfsiidt simidatfons involves predicting multi-dimensional 

and combusting flows. The primary utility or ^ Engineering Sciences, Inc. have 

effects caused by specific injector co gura fcrNASA/MSFC for the past decade. CFD 

been developing a ^Pr o Pri ate > co m P^^ 1 ^^^®^^ f 0 ^^ n ^ I t ^ s ^ period; however , the physica l 
tools and computers have improved dram y J • . in order t0 pro duce useful results, 

submodels used in these analyses must s 1 rem e i emen t s for hydrogen and hydrocarbon fuels, 

Stations of clustery Inland Th.Uc— •> 

which account for real I flu. stroctured t0 repre sent homogeneous and 
codes are based on the FDNS CFD ° model treats the flow as a continuum of 

heterogeneous spray combustion^ T 8 . , thermal or velocity lags between the phases, 

multi-phase, multicomponent fluids w ic vo i ume -of-fluid (VOF) model which represents 

Two heterogeneous models were vaporization, and (2) a Blob 

the liquid core of coaxial or impmger j® h of droplets the size of the injector orifice 

codes to improve- their computational efficiency. 

These spray combustion codes were trsedto — the riuee test ca^which^ the 

subject of the 2nd International Workshop combustion kinetics, 

considered by these investigators to be very '/ahiabfe for code va of hydrogen air 

rocket engine combustion. 

We wish to emphasize that the simulations which we rather 

test the accuracy of the .approximations ^ wefe conducted . Our goal is to 

than represent a definitive analysis ofea ^ h L in rocket engines; hence predicting individual 

accurately predict local temperatures and nmdure F JANNAF standard 

experiments is used for performance prediction with CFD cases, such 

axisymmetric finite-rate (TDK) co p beas to use and of comparable run times for 

"oS^Sctions Secondly, they must provide more detailed predictions of the 
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flowfields near the injector face. Specifically, they must accurately predict the convective mixing of 
injected liquid propellants in terms of the injector element configurations. 


METHODOLOGY 

Homogeneous Spray Combustion Model 

The homogeneous spray combustion CFD codes utilize very general thermodynamics in a 
conventional CFD code. The heterogeneous codes use tabulated properties for the lqui p 
. , . nrnnerties for the vapor phase Thermal and caloric equations of state, vapor pressure, hea 

and transport properties are modeled with the equations of state 

oronosed by Hirshfelder et al 3 ' 4 (we term these the HBMS equations of state) and with conventional 
correlations 5 for the other properties. The property correlations used were not chosen for thei 
absolute accuracy but for their validity over a wide range of temperatures and pressures and for 
requiring a minimum of data to describe a particular species. These correlations are explicit in density 

and temperature. 
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These equations are based on the "theorem of corresponding states" for real fluids, which essentially 
means that the p-v-T relations for all species are similar if these variables are normalized with 
values at the critical point, i.e if reduced values are used. The reduced values in these equations are 
indicated with a subscript r. H. is the ideal gas species enthalpy. Z. is the compresstbil toy ^or a given 
species at the critical point. The HBMS equations are attractive to use because arbitrary correlat o 
for vapor pressure, heat of vaporization, and liquid densities can be used. Since multi-wmponen 
fluid/vapor mixtures may be present in the flowfield, the mixture properties are calculated by the 
additive volume method. This means that multiphase mixtures are treated as ideal solutions 
H2/O2 propellants under conditions where the species become ideal gases, the thermodynamic 

from the CEC code 6 were used. 

The combustion reactions used in the simulations reported herein are shown in Table 1 Not 
all of the reactions were used in all of the combustion simulations. Elementary rate data for these 
reactions are reported by Gardner, et al” Such data are empirical and were obtained for 
hydrogen/air combustion, under conditions far different from those encountered in rocket engines. 
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Table 1. Combustion Model for H 2 /0 2 Reaction 


Chain initiation: 

H 2 + 0 2 = 20H 

1.86H 2 + 0 2 = 1 645 H 2 Q + 0 067 O + 0.142 H + 0.288 OH 

Chain Branching: 

H 2 + OH = H 2 0 + H 
2 OH = H z O + O 
H 2 + O = H + OH 

Q 2 + H = 0 + 0H 

Chain termination: 

O + H + M = OH + M 
2 O + M = 0 2 + M 
2H + M = H 2 + M 

OH + H + M = H 2 Q + M 


The CFD solver used was the Finite-Difference Navier-Stokes code with provision for using 
real fluid properties, the FDNS-RFV code. This code is pressure based; it differs from an ideal gas 
code in the methodology used to relate the pressure correction to the continuity equation and of 
course in the properties subroutines used. The pressure correction (p ; ) equation used in the FDNS- 
RFV code is: 


j ii 

A.( Ui P p p' ) - A • ( p* D P A p') = - A • (p* ui ) - 
p n+, = p n +p' ; P p = y/a 2 ; ui«-D P A P ' 

where the superscripts * and n denote the value at the intermediate and previous time steps, 
respectively. D p is the inverse of the matrix of the coefficients of the convective terms in the finite- 
difference form of the inviscid equations of motion. This is not an obvious definition, but is one 
which has made the FDNS-RFV code a useful solver. The sound speed used in the pressure 
correction equation is that calculated for the real fluid multi-component mixture. 

In all cases simulated, a k-e turbulence model was used to close the mass averaged 
transport equations solved by the code. Our experience is that this incompressible turbulence 
model overestimates the mixing in a combusting flowfield. However, since the liquid propellants 
are also mixed by this model, we concluded that there are currently insufficient data to better tune 
the turbulence model. The homogeneous spray model has been used to simulate: (1) a single 
element like-on-like (LOL) impinger injector element and a single element unlike impinger 
element for the configuration and flow conditions used in the cold-flow experiments; (2) an 
ensemble of injector elements in the Fastrac engine; and (3) several configurations of the vortex 
engine currently being developed. 9 
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Heterogeneous Spray Combustion Model 

Simulations of shear coaxial injector combustion may include models that characterize the 
breakup or atomization of the round liquid jet, subsequent droplet secondary breakup, turbulence 
dispersion, droplet evaporation and gas-phase mixing and combustion. The primary atomization rate 
of the liquid jet is modeled following the work of Reitz and Diwakar . Applications of this model 

to shear coaxial injector test cases, with a volume-of-fluid equation to model the liquid fuel/oxidizer 
jets, were presented by Chen, et al. u . For the present application, since the liquid core length and 
the initial droplet size are specified, the primary atomization model is therefore ignored. 

Particulate Two-Phase Flow Model 

The two-phase interactions are important throughout the life history of the droplets. In the 
initial phase of injection, momentum and energy exchanges through the drag forces and heat transfer 
are dominating. These inter-phase transfer terms appear in the Navier-Stokes equations that are 
solved using the present CFD flow solver. Mass transfer occurs as the particles are heated through 
the surrounding hot gas. Mean gas-phase properties and turbulence eddy properties are used for the 
statistical droplet tracking calculations. 

Droplet Secondary Breakup Model 

The TAB (Taylor Analogy Breakup) model of O’Rouke and Amsden 12 is based on an 
analogy between an oscillating and distorting droplet and a spring-mass system. The restoring force 
of the spring is analogous to the surface tension forces on the droplet surface. The external force on 
the mass is analogous to the gas aerodynamic for ce . The damping forces due to liquid viscosity are 
introduced also based on this model. 

Droplet-Turbulence Interaction 

A two-equation turbulence model is used to characterize the flowfield turbulence quantities, 
such as turbulence fluctuations, eddy life time and length scale. Turbulent effects on particles are 
modeled by assuming the influence of velocity fluctuations on the particles creates statistical 
dispersion of the particles. The velocity fluctuations, which are calculated from the solutions of the 
turbulence kinetic energy, are assumed to follow a Gaussian distribution with standard deviation 
proportional to the square root of turbulence kinetic energy. This magnitude of this statistical particle 
dispersion is then transported following the trajectory of the particles with their radii of influence 
within which coupling effects (also follow the Gaussian distribution) between two phases occur. This 
method is classified as the parcel PDF (cloud) model, by Shang 13 , for turbulent particle dispersion. 
As oppose to the stochastic separated flow (SSF) model, the number of computational particles 
required is drastically reduced for the same statistical representation of the spray. This provides great 
savings in computational effort in performing the spray combustion computations. 

Droplet Evaporation Model 

The droplet evaporation rates and the droplet heat-up rates are determined using the general 
evaporation model of Schuman 14 , which is continuously valid from subcritical to supercritical 
conditions. This vaporization model was extended from the classical approach 1S , by neglecting the 
effects of solubility of the surrounding gas into the droplet. However, this approach satisfies the 
global transient film continuity equation for the drop vapor and the ambient gas to obtain the 
expressions consistent for the molar flow rates, 
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Chemical Reaction Model 

A finite-rate chemistry model with point-implicit integration method is employed in the 
present study. A 9-reaction kinetics model of Anon 16 is used for modeling the H 2 -O 2 combustion. 
The initiation reaction used produced OH. This chemistry model is listed in Table 1. 


SIMULATIONS OF THE RCM-1 EXPERIMENTS 

The LN2 cases, RCM-1 -A and -B, were simulated with the homogeneous spray model. The 
flow predicted resembles a dense fluid jet with strong density gradients in the shear layer. Such a 
flow has been observed in a similar super-crit ical n itrogen jet experiment reported by Chehroudi, et 
al 17 . These predictions should compare well to the DLR experimental data. If the comparisons are 
not good, adjustment of the parameters in the two-equation k-e or the initial turbulence level 
parameters could be made for a better fit of the data . Such tuning has not previously been made since 
appropriate test data were not available. For a definitive analysis of the experiments, conjugate heat 
transfer to the injector hardware and consideration of the duration of the experiment should be made. 
The jet is discharging into a gaseous nitrogen environment, the recirculated gas should become slowly 
cooled until a steady state is reached. Since the temporal variation of the recirculating gas 
temperature was not reported, the time that the CFD simulation should be terminated can not be 
determined. Since the measurements were made very close to the injector exit, good simulation of 
the gas temperature might not be crucially important. 

The injector configuration and flow conditions for the cryogenic nitrogen jet of the RCM-1 
test cases are illustrated in Fig. 1. It can be seen that the chamber pressure for both cases is above 
the critical pressure of nitrogen. A 101x1 1-mesh system was used to discretize the injector section, 
while the chamber section was modeled by a 301x10 1-mesh system for Case RCM-1 -A.. The same 
grid system was used to simulate both RCM-l-A and RCM-1 -B test cases. The numerical result of 
RCM-1 -A test Case at the locations specified by IWRCM was plotted as shown in Figures 2-6. 
Notice the temperature profiles in Figure 4. These two cold flow cases are not steady-state, although 
the simulations assumed this to be the situation. The simulations presented represent a time-slice at 
some arbitrary time. Figure 7 shows the flowfi eld near the injector tip. A finer gnd system (101x15, 
and 301x141) was employed to simulate the RCM-1 -B. The numerical results of RCM-l-B test case 
are plotted in Figures 8-12. The flowfield is presented in Fig. 13. Notice that only a small segment 
of the chamber is shown so that gradients in the flowfield may be clearly seen. 
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(a) Streamwise Velocity 


(b) Density 


(c) Turbulence Kinetic Energy 




Figure 2 Flow Properties at the Injector Exit of RCM-l-A 
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Figure 4 Temperature Profiles at Various Streamwise Locations of RCM-l-A 
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Figure 6 Turbulent Kinetic Energy Profiles at Various Streamwise Locations of RCM-l-A 
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Figure 7 Flow Properties Near the Injector of RCM-l-A 
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Figure 8 Flow Properties at the Injector Exit of RCM-l-B 
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Figure 9 Density Profiles at Various Streamwise Locations of RCM-l-B 
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Figure 12 Turbulent Kinetic Energy Profiles at Various Streamwise Locations of RCM-l-B 
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Figure 1 3 Flow Properties Near the Injector of RCM- 1 -B 
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SIMULATIONS OF THE RCM-2 EXPERIMENTS 

The sub-critical combustion case, RCM-2, was simulated with both the heterogeneous and 
the homogeneous spray combustion models. The MASCOTTE test data should be better than any 
which have been previously used to tune the several parameters in these models. It is unreasonable 
to expect that spray flames, even of hydrogen and oxygen, can be accurately predicted without 
extensive model validation with test data representative of the conditions which exists in rocket 
engine combustion chambers. Even global data like chamber pressure and thrust have not been 
obtained for single coaxial element combustor flows. The IWRCM data provide a good starting 
point, but no CFD model tuning has yet been attempted for such experiments. Direct comparisons 
of predictions to test data at this point will not establish which of several modeling techniques is best. 

The MASCOTTE single injector test chamber was used in a series of experimental programs 
for subcritical and/or supercritical H 2 -O 2 combustion. In the subcritical spray combustion test case 
(RCM-2), the designed chamber pressure is 10 bar (or 9.87 atm). The injector orifice diameter for 
the liquid oxygen (LOX) injection is 5 mm surrounded by an annular gaseous hydrogen jet with 
channel width of 6.4 mm. The overall O/F ratio for this case is 2. 1 1 (see the test conditions given in 
Table 2). 


Table 2. RCM-2 Test Case Operating Conditions 


Conditions 

h 2 

0 2 

Pressure 

1 MPa 

1 MPa 

Mass flow rate 

23.7 g/s 

50 g/s 

Temperature 

287 K 

85 K 

Density 

0.84 kg/m3 

1170 kg/m 3 

Cp 

14300 J/kg/K 

1690 J/kg/K 

Velocity 

319 m/s 

2.18 m/s 

Viscosity 

8.6E-6 kg/m/s 

1.94E-4 kg/m/s 

Surface Tension 

- 

1.44E-2 N/m 


The computational model includes the injector geometry, the combustion chamber and the 
nozzle section. A 10-block structured mesh is generated (the total number of grid points equals 
14,444) for the two-phase flow computation. Relative high grid density (about 10 micron spacing)is 
packed in the injector lip region for the purpose of better flow resolution and flame holding in the 
expected area. The LOX core length of 7.8 mm is assumed, which serves as the particle injection 
boundary with the fixed particle size (82 microns), velocity (10 m/s) and angle distributions given in 
the problem specification. Fixed mass-flow boundary conditions are used at the inlet while all flow 
properties are extrapolated at the nozzle exit. Supersonic exit flow develops as part of the solution. 

The computation starts with a cold flow with inlet and chamber pressure specified. The two- 
phase flow particle breakup and evaporation model models are activated from the beginning. The 
time step size of the time-marching solution method is 1 psec. After 1000 time steps of cold flow 
run, a heat source is introduced in the lip region between and oxygen and hydrogen streams where 
a recirculation zone is established. At the same time, the finite-rate chemistry model is turned on to 
start the flame spreading throughout the chamber. The chamber pressure drops at the beginning until 
the flame fills up the entire chamber. Then, the pressure started to build up to the expected level 
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when the inlet and exit flows show satisfactory mass conservation condition. The calculated averaged 
chamber pressure is around 9.96 atm. The majority of the LOX particles do not survive very far 
downstream of the injector exit. Some particles along the chamber axis do survive up to 70 mm 
downstream of the injector. 

The time-averaged temperature, temperature standard deviation, species mass-fraction 
contours and temperature profiles at specified locations are plotted in the following figures. These 
data are prepared as requested for data comparison purpose. 

Figure 14 shows the mean temperature and standard deviation through the entire length of 
the combustion chamber. A close up view of the nozzle tip region is also shown in this figure. Figure 
15 shows the OH and O 2 and Figure 16 the H 2 and H 2 0 concentration profiles, respectively, in this 
same region. Figures 17-25 show radial temperature and standard deviation profiles at various axial 
locations. Figures 26-30 show the axial temperature profiles at various radial locations. Figure 30 
shows this profile at the near wall location. The flame predicted with this model is long and narrow. 
The recirculation zone is very long. 

The RCM-2 experiment was also simulated with the homogeneous spray combustion model. 
The volume upstream of the injector element tip was neglected for this simulation. The grid use for 
the internal element flow was 61 X 43; for the chamber it was 301 X 101. The nozzle was not 
simulated. This grid system had a minimum grid spacing of 60 microns in the wake behind the lip 
separating the LOX and hydrogen streams. The boundary conditions used are shown in Figure 3 1 . 
An equilibrium and several finite rate solutions were obtained for this configuration. The rate of the 
global initiation reaction was set fast enough to stabilize the flame near the start of the shear layer. 
This rate also essentially eliminated the waviness in the shear layer separating the LOX and hydrogen 
streams, without averaging the solution. The stoichiometric coefficients in the global rate expression 
were determined by an equilibrium calculation for a stoichiometric flame at the expected chamber 
pressure. Such a practice produces temperatures with one rate expression which are very close to 
those resulting from using a more detailed reaction mechanism. 

The equilibrium solution at the interface between the internal element flow and the flow at the 
nozzle tip are shown in Figure 32. The temperature profiles in the radial and axial directions are 
shown in Figures 33 and 34, respectively. The temperature and oxygen and OH concentration 
profiles are shown in Figure 35. The wall temperature profile is shown in Figure 36. All of these 
figures are for the equilibrium solution. The finite rate solutions for the single global reaction and for 
the global plus the elementary reactions of Table 1 were also obtained. As expected, the finite rate 
solutions were slightly cooler than the equilibrium solutions. The predictions are very similar in all 
results to those just shown. To emphasize this point the wall temperature profiles for all three cases 
are shown in Figure 36. Even though the global rate was set fast enough to stabilize the flame with 
this grid system, it was not so fast that equilibrium conditions were obtained. 

Comparing the heterogeneous and homogeneous solutions, the former produced a longer, 
thinner flame than the latter. Parameters in the spray combustion model could have been set such that 
the solutions matched very closely, or so that both could match test data. Such a step cannot be made 
until the RCM test data are published and the CED models tuned. An optimum rocket engine spray 
combustion model cannot be determined until this next validation step is undertaken. 
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(c) Time- Averaged Temperature (K) for X up to 150 mm 


Figure 14 Time- Averaged Temperature and Temperature Standard Deviation of RCM-2 
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(a) Time- Averaged OH Mass-Fraction Contours for X up to 1 50 mm 




(b) Time- Averaged 02 Mass-Fraction Contours for X up to 150 mm 


Figure 15 Time- Averaged OH and O 2 Mass Fractions of RCM-2 
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Figure 16 Time- Averaged FT 2 and H 2 0 Mass Fractions of RCM-2 
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Figure 17 Radial Profiles of Mean Temperature and Standard Deviation at X/D = 2 of RCM-2 
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Figure 18 Radial Profiles of Mean Temperature and Standard Deviation at X/D =10 of RCM-2 
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Figure 23 Radial Profiles of Mean Temperature and Standard Deviation at X/D = 43 of RCM-2 
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Figure 24 Radial Profiles of Mean Temperature and Standard Deviation at X/D = 50 of RCM-2 
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Figure 25 Radial Profiles of Mean Temperature and Standard Deviation at X/D = 60 of RCM-2 
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Figure 26 Axial Profiles of the Mean Temperature at Y/D = 1 of RCM-2 
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Figure 27 Axial Profiles of the Mean Temperature at Y/D = 2 of RCM-2 
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Figure 28 Axial Profiles of the Mean Temperature at Y/D = 3 of RCM-2 
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Figure 30 Axial Profiles of the Mean Temperature Near the Outer Wall of RCM-2 
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Figure 34 Axial Profiles of the Mean Temperature at Various Radial Locations of RCM-2 
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Figure 35 Temperature and Species Concentrations Near the Injector of RCM-2 
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Figure 36 Near Wall Temperature Distributions for Various Chemistry Models of RCM-2 



SIMULATIONS OF THE RCM-3 EXPERIMENT 

The super-critical combustion case, RCM-3, was simulated with the homogeneous spray 
combustion model. Any drops present will be highly unstable; therefore, this model should represent 
the flow rather well. Local equilibrium and simplified finite-rate combustion submodels were used 
and the results for the two simulations compared well. More detailed combustion submodels were 
attempted, but proved to behave too poorly for successful simulations. 

The preponderance of super-critical spray combustion models which have been reported have 
been extensions of sub-critical models. Such models encounter a basic problem in over emphasizing 
the role of surface tension. Since surface tension is zero for super-critical conditions, drops should 
not exist. Although such drops can be observed experimentally, they are extremely unstable and do 
not survive very long. The homogeneous CFD model was developed to account for the major 
physical effects which do exist. Namely, the large density and momentum differences which exist in 
multi-phase super-critical flows. Such a model allows one to accurately relate the inlet conditions 
at the injector face to boundary conditions for the CFD simulation. This relationship is essential to 
predicting the effects of injector element configuration and inlet momentum vector on the convective 
mixing and cross winds which occur in practical rocket engines. Otherwise, one is forced to use the 
historical method of creating costly experimental data bases from which to choose designs. 

The injector configuration and flow conditions for the supercritical combustion of the 
RCM-3 test case are presented in Fig. 37. This is uni-element shear coaxial injector with LOX 
and GH 2 propellants. The numerical simulati on wa s conducted with some simplification because, 
initially, detailed information was unavailable; such as: (1) the flare of LOX injector near the exit 
was neglected; (2) the injector was flush at the chamber head-end instead of protruding into the 
chamber because the outer diameter of hydrogen tube and distance between the chamber head- 
end and the injector exit were not known; (3) the nozzle was not included because of insufficient 
information about the chamber tail-end and nozzle geometry; and (4) the coolant (later found to 
be helium) for the chamber wall was not included because its flow rate and properties were not 
specified. As can been seen, the chamber pressure (60 bar) is well above the critical pressure of 
oxygen; hence, the homogeneous real-fluid model was used to simulate this test case. A two- 
zone mesh system (61x39 and 301x101) was used to model the injector section and the 
combustion chamber. 


The combustion reactions in this high pressure experiment are expected to be in local 
thermodynamic equilibrium and were simulated as such. To demonstrate the methodology, two 
finite-rate simulations were also made with a subset of the reactions in Table 1. The single global 
reaction which produces radicals as well as water provides a good estimate of the temperature field. 
Its rate was set to attach the flame near the Injector tip. Since the radicals are not rigorously 
simulated with the single reaction, a second finite-rate simulation was made with the 2-body reactions. 
Backward reaction rates are determined with eq uilibrium constants. For high pressure cases such 
combustion modeling is essential to keep the computation stable. 

The chemistry and turbulence models used in our simulations do not make use of probability 
density functions (PDFs) because most of the shear layers formed by the injector element should be 
continuum. The only regions for which this might not be the case are the intermittent edges of the 
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shear layers. Pope 8 terms these regions the "viscous superlayer". The thickness of these layers are 
inversely proportional to the Reynolds number to the 0.75 power. For these high speed coaxial jets, 
they should be very thin. 

The flow predicted at the injector tip is shown in Figure 38. The radial temperature profiles 
predicted at several axial stations are shown in Figure 39. The axial profiles at several radial locations 
are shown in Figure 40. The temperature and o xy gen and OH concentration profile fields are shown 
in Figure 41. The combustion models used do not predict chemiluminescent OH, which might be 
observed in the experiments. These results are shown for the equilibrium combustion model. Results 
for the two finite-rate combustion simulation s are very similar, hence they are not shown. The wall 
temperature distributions for all three cases are compared in Figure 42, and as noted the results are 
very similar. 



Figure 37 Configuration of the RCM-3 Case (Homogeneous Spray Model) 


(a) Axial Velocity 



e 

E 


8 

S 

00 

a 


■o 


£ 


(b) Density (c) Temperature 



Figure 38 Flow Properties at the Injector Exit of RCM-3 (Homogeneous Spray Model) 
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Figure 39 Radial Profiles of Mean Temperature at Various Axial Locations of RCM-3 



Figure 40 Axial Profiles of Mean Temperature at Various Radial Locations of RCM-3 
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Figure 41 Temperature and Species Concentrations Near the Injector of RCM-3 
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Figure 42 Near Wall Temperature Distributions for Various Chemistry Model of RCM-3 
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CONCLUSIONS 


The following conclusions were drawn from performing CFD simulations of the three RCM 
test cases for the 2nd IWRCM. 

1 . A homogeneous and a heterogeneous spray combustion CFD model have bee developed to 
simulate combustion in rocket engines. Since neither of these models is expected to be accurate 
until critical parameters are evaluated from test data, simulation comparisons to the MASCOTTE 
type experiments are needed. 

2. The utility of either CFD model cannot be determined until values of critical parameters are 
determined and efforts to optimize the computational efficiency of the models are performed. 


3. Although the CFD rocket engine models provide much more detailed information concerning the 
vaporization, mixing, and combustion process, their place in the design process is yet to be 
identified. Older more approximate rocket "performance" models are difficult to displace. 
Furthermore, every physical process thought to be present in the engine does not have to be 
modeled to create a useful design code. There are more knobs to adjust in the code than there 
are experimental data to justify their turning. 

4. The experiments conducted in preparation for the 2nd IWRCM appear to be a significant first step 
in providing test data valuable to CFD modelers. However, blind comparisons of CFD model 
predictions to such data are premature. The CFD modelers have not previously had sufficient test 
data properly specify the many assumptions which are necessary to simulate such complex flows. 

5. Better communication between analysts and experimenters needs to be accomplished. Can the 
modeler simulate the experiments which are being performed? Can the data obtained from the 
experiment critically test the model? 
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